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ABSTRACT 

Context. Swift observations suggest that the X-ray afterglow emission of some gamma-ray bursts (GRB) may have 
internal origins, and the conventional external shock (ES) cannot be the exclusive source of the afterglow emission. 
Aims. If the central compact objects of some GRBs are millisecond magentars, the magnetar winds could play an 
important role in the (internal) X-ray afterglow emission, which is our focus here. 

Methods. The dynamics and the synchrotron radiation of the termination shock (TS) of the magmnetar winds, as well 
as the simultaneous GRB ES, are investigated by considering the magnetization of the winds. 

Results. As a result of the competition between the emission of the wind TS and the GRB ES, two basic types of X-ray 
afterglows are predicted, i.e., the TS-dominated and the ES-dominated types. Moreover, our results also show that both 
of the two types of afterglows have a shallow-decay phase and a normal-decay one, as observed by the Swift satellite. 
This indicates that some observed X-ray afterglows could be (internally) produced by the magnetar winds, but not 
necessarily GRB ESs. 
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1. Introduction 

The conventional fireball model for gamma-ray bursts 
(GRBs; see reviews by Piran 2005; Meszaros 2006) sug- 
gests that the burst emission is produced by internal dis- 
sipations in a relativistic, expanding fireball, while the af- 
terglow emission arises from the deceleration of the fireball 
in the circumburst medium by exciting an external shock 
(ES). However, Swift observations have showed that some 
X-ray afterglow emission components may sometimes also 
have internal origins. The most solid evidence of this is 
' the rapid rise and decline of X-ray afterglow flares (e.g., 
Burrows et al. 2005; Yu & Dai 2009). Meanwhile, the ES 
model may also be challenged by some plateau emission 
that is followed by an abrupt cutoff (Troja et al. 2007). 
These phenomena indicate that internal afterglow emis- 
siorQ could be more usual than considered before, at least 
in the X-ray band, whereas the ES is only one choice among 
various afterglow origins. In this case, some observed chro- 
matic breaks in the afterglow light curves may also become 
more understandable. The long-lasting internal emission 
strongly requires long-lasting activities (energy release) of 
GRB central compact objects after the bursts, and to a 
certain extent, spinning-down magnetars are very plausible 
candidates (e.g., Dai et al. 2006; Troja et al. 2007). 



1 In this paper we use the word "internal afterglow emission' 
to generally refer to the emission from any sites inside the ES. 



It has been proposed for nearly two decades that GRBs 
could be produced by spinning- down, millisecond magne- 
tars (Usov 1992; Paczynski 1992; Duncan & Thompson 
1992). In this case, the relativistic fireball responsible for 
a GRB could be associated with a neutrino-driven wind 
(Thompson 1994; Metzger et al. 2007; Bucciantini et al. 
2007, 2009), with a magnetic reconnection-accelerated wind 
(Drenkhahn 2002; Drenkhahn & Spruit 2002; Giannois 
2008), or with hyperaccretion onto the magnetars (Zhang 
& Dai 2008, 2009). Of more interest here, the postburst 
magnetars could still keep abundant rotational energy, al- 
though a considerable fraction of the energy has been ex- 
pended on the bursts. With the spin-down of the magne- 
tars, this remaining energy can be released and drive a con- 
tinuous energy outflow, i.e., a postburst wind. In view of the 
huge energy carried by the wind, Yu et al. (2010) suggested 
that some observed X-ray afterglows could be directly con- 
tributed by such winds, according to a qualitative compar- 
ison between the model and the observations. However, in 
that work, no specific radiation mechanism was addressed 
for the wind emission besides introducing a free constant 
for the radiation efficiency of the winds. 

In this paper we try to search for a possible specific 
mechanism for the wind emission. Following Dai (2004) and 
Yu & Dai (2007), we suggest that the termination shock 
(TS) of a magnetar wind, which arises from the interaction 
of the wind with the GRB fireball and the surrounding 
medium, may play an important role in the wind emission. 



2 



Z. Mao et al.: GRB X-ray afterglows from magnetar winds 



Based on a more careful consideration, we redescribe the 
emission from the wind TS and the GRB ES by considering 
(1) the existence of the TS after the spin-down timescale 
and (2) the magnetization of the wind. 

2. The spin-down of a magnetar and the magnetar 
wind 

Some short-term and violent processes mean that a newly 
born magnetar could drive an energetic, relativistic, com- 
plicated outflow to produce a GRB. During the first tens 
of seconds after the birth of the magnetar, its spin evolu- 
tion must therefore be very complicated. In contrast, on 
the relatively longer (afterglow) timescales concerned here, 
the spin-down of the magnetar would be mainly controlled 
by the electromagnetic torque (leading to an MHD outflow) 
and the torque connected with gravitational wave radiation 
(Zhang & Meszaros 2001). For the latter, if the equatorial 
ellipticity of the magnetar is mainly determined by the stel- 
lar magnetic field, the energy losses through gravitational- 
quadrupole radiation would not be essential for the spin- 
down of the star unless the mag netic field is > 3 x 10 17 G 
(Usov 1992). 

Therefore, for a magnetar with a magnetic field of 
~ 10 14 ~ 15 G, its postburst spin-down luminosity can be 
simply estimated by (Shapiro & Teukolsky 1983) 
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where R sfi = i? s /10 6 cm, B sM = B s /10 14 G, and /'„.,. , = 
P s ,i/1 m s are the radius, magnetic field strength, and initial 
spin period of the magnetar, respectively. The characteristic 
spin-down timescale reads as 
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where J S) 4g = 7 S /10 45 g cm is the stellar moment of iner- 
tial. The time t is defined in the observer's frame, but the 
time zero is actually not strictly set at the GRB trigger 
time since Eq. (fTJ) is only valid about 10 — 100 s after the 
birth of the magnetar. In other words, t = in this pa- 
per approximately represents the transition time from the 
prompt phase to the afterglow phase. Such a shift in the 
time zero would not influence the calculation of the after- 
glow emission on much longer timescales (t > 10 2 ~ 3 s). 

With the release of the stellar rotational energy after 
the burst, an MHD outflow is produced whose luminos- 
ity is determined by the spin-down luminosity as shown in 
Eq. ([1]). Unlike the GRB fireball that could be a neutrino- 
driven wind, this postburst outflow is likely to be initially 
Poynting flux-dominated, mixed by an electron-positron 
plasma. Since the fluctuating component of the magnetic 
fields in the outflow can in principle be dissipated by 
magnetic reconnection, the associated leptons are acceler- 
ated and the outflow gradually becomes an ultrarelativistic 
kinetic-energy flow consisting mainly of electron-positron 
plasma (lepton-dominated wind; Coroniti 1990). For the 
bulk Lorentz factor of the accelerated wind, we here arbi- 
trarily take T w ~ 10 4 as a fiducial value following Atoyan 



(1999), who argues that the Crab pulsar wind had Lorentz 
factor T w ~ 10 4-7 to interpret the measured radio spec- 
trum of the Crab Nebula. 



3. The dynamics 

Owing to its ultrarelativistic velocity, the magnetar wind 
catches up to and collide with the decelerating GRB ejecta. 
The interaction of the wind with the ejecta and the sur- 
rounding medium gives rise to a relativistic wind bubble, 
which includes four different regions separated by a con- 
tact discontinuity surface and two shocks (i.e., the TS of 
the wind that propagates into the cold wind and the GRB 
ES that propagates into the surrounding medium^). To be 
specific, the four regions are (1) unshocked medium, (2) 
shocked medium including the GRB ejecta, (3) shocked 
wind, and (4) unshocked cold wind close to the front of 
the TS. The shocked regions are emission regions. 

We denote some quantities of region i as follows: n[ is 
the particle (proton or electron) number density, B[ the 
magnetic field strength, e[ the internal energy density, P! 
the pressure, where the primes refer to comoving frame, and 
the subscript i represents regions 1—4. I\, = (1— T^ 2 ) 1 / 2 , 
and Ui — PiTi are the bulk Lorentz factor, velocity, and 
four- velocity, respectively, measured in the local medium's 
rest frame. It is convenient to define a parameter 



B , 



Ann' 4m e c 2 



to denote the degree of magnetization of the unshocked 
wind near the front of the TS, where m e is the electron 
rest mass and c the speed of light. By assuming the wind is 
isotropic, the comoving electron density of the unshocked 
wind can be expressed as n' 4 = L m( j / '[4ir R 2 r 2 m e c 3 (1 + a)} 
with L4 = T w , and R is the radius of the bubble in the thin 
shell approximation^. As considered in Zhang & Kobayashi 
(2004) and Fan et al. (2004), a possible slow decrease in B'\ 
at large radii due to magnetic reconnection can be ignored. 
Thus same as n.4, B'\ would approximatively decrease as 
R~ 2 with the expansion of the bubble, and thus the param- 
eter a can be regarded as a constant. 

To connect the properties of the two sides of the shocks, 
some shock jump conditions need to be listed as (Kennel & 
Coroniti 1984; Zhang & Kobayashi 2004) 



r. 



= 4L 2 
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(4) 
(5) 
(6) 
(7) 



2 The initially-formed reverse and forward shocks in the GRB 
ejecta would cross the ejecta very quickly (~ 100 s), so we do 
not consider the effect of these two shocks. 

3 Following Eq. (7) in Dai (2004), the thickness of the relativis- 

x 



tic wind bubble, A = R 



(> + ^)(> 



- 1 



< R, 



could be negligible with respect to the bubble radius, where the 
similarity variable \ i s defined in Eq. (|15|l . 
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where an adiabatic index of 4/3 is adopted, m p is the proton 
rest mass, T 34 = (T 3 /T 4 + T 4 /T 3 ) /2 w r 4 /(2r 3 ) is the 
bulk Lorentz factor of region 3 measured in the rest frame 
of region 4, 
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and u 3 t is the four- velocity of region 3 measured in the rest 
frame of the TS. When T 34 >• 1, which is easily satisfied 
in this model (Yu & Dai 2007), Eqs. (9) and (10) can be 
simplified as 
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and u 3 t can be solved analytically from 8(a+l)u 3t — (8a 2 + 
10cr + l)u\ t + a 2 = (Zhang & Kobayashi 2004). 

Following Eq. (4), the total kinetic energy of region 2 
can be calculated as Ek, 2 = (r2 — l)(m j+m sw )c 2 + r 2 (r 2 — 
l)m sw c , where m e j and the rest masses of the GRB 

ejecta and swept-up medium, respectively. Energy conser- 
vation requires that any increase in Ek.2 should be equal to 
the work done by region 3, i.e., 



dE k ,2 =SW = 4nR 2 (P; K3 + P Bj3 )dR, 



(13) 



where the thermal and magnetic components of the 
pressure of region 3 can be respectively calculated by 



PL * = = \T\^m e (?f a f h and P' B 3 = S' 2 /(8tt) 



th,3 — 3 C 3 — 3 J 

8r 34 n' 4 m e c 2 af 2 according to the shock jump conditions. 
Substituting the expression of Ef.,2 into Eq. (13), we can 
obtain 
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Meanwhile, the relationships between the properties of 
the two sides of the contact discontinuity surface show 
(Blandford & McKee 1976) 
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where the similarity variable can be solved to be 
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For a closing dynamic equation set, we also introduce an- 
other two equations in order to calculate the increasing 
masses of regions 2 and 3 as 



dm s 



dR 



= A / nR 2 nim p 



(18) 



and (Dai & Lu 2002) 



dm 3 
~dR 



AttR z 



r 34 /? 34 
r 3 /? 3 



n 4 m e . 
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Finally, we would like to pointed out that, in Dai (2004) 
and Yu & Dai (2007), the authors only took the existence 
of the TS before T m into account since the shallow decay 
afterglow phase during ~ 10 3-5 s was mainly concerned 
there. In contrast, here such an artificial cutoff of the TS 
has been abandoned for a more general investigation. 

Combining the above dynamic equations with the re- 
lationships of dt = (1 — (5 3 )dR/ (J3 3 c) for region 3 and 
dt = (1 — /?2)d-R/(/3 2 c) for region 2, we can numerically 
calculate the temporal evolution of the shock dynamics. As 
mentioned in Sect. 2, since t = is set at the transition 
time from the prompt phase to the afterglow phase, the 
initial conditions for our calculation is taken as Ri = 10 16 
cm, T2,i = 150, and £7^2,1 = 10 51 erg, which correspond 
to the deceleration timescale of the GRB ejecta. Figure 1 
shows the evolution of the bulk Lorentz factors of the wind 
TS and the GRB ES with different values of a. The nearly 
overlapping of the curves with a varying from to 3 indi- 
cates that the dynamic evolutions of both the shocks are 
insensitive to the degree of magnetization of the wind. It is 
not surprising to obtain such a result, which actually has 
been know for a long time for ordinary pulsar wind nebulae 
(e.g., Emmering & Chevalier 1987; Bucciantini et al. 2003; 
Del Zanna et al. 2004). 

4. The radiation from the shocks 

With the propagation of the shocks, the bulk kinetic en- 
ergy of the GRB ejecta and the magnetar wind would be 
gradually converted into internal energy of the shocked 
materials. As usual, we assume that the internal energy 
of the shocked medium is shared by magnetic fields, elec- 
trons, and protons with fractions es, e e ~ ^JtB (Medvedev 
2006) and 1 — e e — e B , respectively. For the shocked wind, 
the energy fractions of the leptons and magnetic fields can 
be determined by the shock jump condition with a cer- 
tain a. Considering the shock acceleration of charge parti- 
cles, the electrons in shocked regions are assumed to dis- 
tribute as n , , cx with a minimum electron Lorentz 



(14) factor: 



7. 



e g p (m p /m e )(T2 — 1) for region 2 and 



imfl = g P ey(n' 3 m e c 2 ) for region 3 with g p = (p-2)/ 
Owing to the synchrotron cooling of the electrons, a cool- 
ing Lorentz factor j' ci = 6irm e c/(a T B' 2 Tit) should also 
be defined by equaling the cooling time to the dynamic 
time, where a T is Thomson cross section. Characterized 
by j' m , j' c and p (the subscript i is omitted hereafter), a 
quasi-static distribution of the electrons can be written as 
(Huang et al. 2000) 



ny oc 



7'l < i < 1'h, 

where 7^ = max^, 7^), 7^ = min^^), and x = 2 for 
l'c < 7m and x = P for ic > 7m- 7m ~ q e B'R/ (m e c 2 ) is 
the maximum electron Lorentz factor with q e the electron 
charge. 

For the X-ray band of interest in this paper, we only 
consider the synchrotron radiation of the electrons in both 
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regions 2 and 3. The synchrotron emission coefficient at 
frequency v' can be given by (Rybicki & Lightman 1979) 



from the numerical results shown in Fig. 1. Combining with 



where 



1 

47T 



7 A-J 



7' 



Ti, 



VSq 3 e B' 



Ki{y)dy 



di, (21) 



3g e 5'7 ,2 /(47rTO e c), Ks(y) is the Besscl func- 
tion. Then by assuming the viewing angle is zero, we can 
calculate the synchrotron flux density received by the ob- 
servers as (Yu et al. 2007) 



F v {t) 



1 

Am 



sin 9 



T 3 (l~/3 cos 6») 3 2 



cos 9d9, 



(22) 



where d^ <~ 10 28 cm is the luminosity distance of the 
GRB and V' is the comoving volume of the emission re- 
gion. Considering the equal-arrival-timc effect for a fixed 
observer's time t, different angles 9 should correspond to 
different radii R, satisfying (Huang et al. 2000) 



t 



f R L 

Jo 



/3(r) cos ( 



/9(r)c 



■dr. 



(23) 



Therefore, as functions of radius, j' v , and V should also be 
functions of 9 in the integral in Eq.(22). 

Some X-ray (^1 keV) light curves contributed by the 
wind TS and the GRB ES as well as the combination of 
them are shown in Fig. 2 with different values of a and e e . 
These two parameters control the energy densities of the 
electrons and magnetic fields, and thus determine the ra- 
diation efficiency of the two emission regions. Fig. 2 shows 
that the emission from the ES is mainly dependent on the 
value of e e (upper panel) but insensitive to a (lower panel). 
This is because the variation in a has little influence on 
the shock dynamics. In contrast, the lower panel of Fig. 
2 shows that the emission from the TS is strongly depen- 
dent on the value of a. More specifically, because a weak 
magnetic field can suppress the synchrotron radiation and 
a strong magnetic field reduces the energy fraction of the 
electrons, neither too low nor too high a leads to strong TS 
emission, which can be seen clearly in Fig. 3. For the model 
parameters adopted in Fig. 3, a ~ a* 0.05 results in the 
strongest TS emission. When the value of a is much closer 
to a* or e e is sufficiently small (e.g., e e <~ 0.02), the TS 
emission is very easy to exceed the ES emission. As a re- 
sult, three types of X-ray afterglows can be predicted, i.e., 
the TS-dominated type, the ES-dominated type and the 
intermediate type, as shown by the left, right, and middle 
panels of Fig. 2, respectively. Qualitatively, such a result is 
basically consistent with the conjecture proposed by Yu et 
al. (2010). 

As found previously (e.g., Dai 1998a, b; Zhang & 
Meszaros 2001), the profile of the ES-dominated afterglow 
light curves (e.g., the cases in the right panel of Fig. 2) is 
in good agreement with the observational light curves in- 
cluding a shallow decay phase and a normal decay phase 
(Zhang et al. 2006), because of the energy injection into 
the ES before T m . More surprisingly, on the other hand, 
we find that a very similar profile also appears in the TS- 
dominated cases (e.g., the cases in the left panel of Fig. 2). 
Please note that the "normal decay" here is actually not 
produced by the ES. To understand such an unexpected 
"normal decay" , we adopt an approximative dynamic evolu- 
tion of the TS as T 3 oc i -1 / 2 for t > T m , which is extracted 



'1/2 



Tm,3 



OC t 



1/2 



oc L m d/(r|t) oc t , we can get B' 3 oc t 
and Y c 3 oc i 1 / 2 . Then following Sari et al. (1998), we can 
easily derive F v oc i -1 from the characteristic spectral 
quantities v m oc t°, v c oc t° , -F^, max oc and the rela- 
tionship v m < ux < v c . The ratio of F v /L m ^ oc t indi- 
cates that the X-ray luminosity of the wind emission docs 
not track the spin-down luminosity of the magnetar with 
a constant fraction, which is assumed in Yu et al. (2010). 
However, the more important conjecture proposed by Yu et 
al. (2010) that some observed X-ray afterglows are directly 
produced by the magnetar winds can be still strongly fa- 
vored by the good agreement between the profiles of the 
TS-dominated afterglow light curves and the observational 
ones. 

Finally, although both the ES- and TS-dominated after- 
glow light curves are qualitatively consistent with the ob- 
servations, some tiny differences between them could also 
be found by some careful statistical studies. For example, 
the shallow decay phase of the TS-dominated afterglows 
would be slightly flatter than for the ES-dominated after- 
glows, and the same is true for the so-called normal decay 
phase. 



5. Summary 

Following Dai (2004) and Yu & Dai (2007), we investigated 
the emission from a TS of a magnetar wind both before 
and after T m as well as the simultaneous emission from the 
GRB ES, where the magnetization of the wind was taken 
into account. The competition between the emission of the 
two shocks is mainly determined by the model parameters 
a and e e . As a result, an ES-dominated type and a TS- 
dominated type of afterglows are predicted, both of which 
are roughly consistent with the observational X-ray after- 
glows. This indicates that some observed afterglows could 
come from the wind TS, but not necessarily from the ES. 
To a certain extent, we regard such TS emission as general 
internal afterglow emission, although the radius of the TS 
is nearly as large as that of the ES. To distinguish between 
these two types of afterglows, some tiny differences in the 
temporal indices could be useful for some statistical studies 
such as that in Yu et al. (2010). In this paper, very high val- 
ues of a were not considered, because in those cases the TS 
may disappear (Giannios et al. 2008) and the injecting flow 
is probably dominated by Poynting flux. This was studied 
thoroughly in Dai & Lu(1998a,b) and Zhang & Meszaros 
(2001). 
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Fig. 3. The cr-dependence of the synchrotron flux of the 
GRB ES (dashed line) and the wind TS (solid line) at the 
time of 10 4 s. The vertical dotted line represents a* — 0.05. 
The model parameters are the same as those in Fig. 2. 



